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Abstract 

The catalysis of chiral symmetry breaking by a magnetic field in the massless weak- 
coupling phase of QED is studied. The dynamical mass of a fermion (energy gap in the 
fermion spectrum) is shown to depend essentially nonanalytically on the renormalized 
coupling constant a in a strong magnetic field. The temperature of the symmetry 
restoration is calculated analytically as Tc ~ mdyn, where nidyn is the dynamical mass 
of a fermion at zero temperature. 
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In this talk I will discuss dynamical chiral symmetry breaking in a magnetic field and at 
finite temperature. The talk is based on the recent papers with V.Miransky and I.Shovkovy 

The dynamics of fermions in an external constant magnetic field in QED was considered 
by Schwinger long ago 0]. In that work, while the interaction with the external field was 
considered in all orders in the coupling constant, the quantum dynamics was treated pertur- 
batively. There is no dynamical chiral symmetry breaking in QED in this approximation. 
Also, chiral symmetry breaking is not manifested in the weak coupling phase of QED in the 
absense of a magnetic field, even if it is treated nonperturbatively. We will show that a con- 
stant magnetic field B changes the situation drastically, namely, it leads to dynamical chiral 
symmetry breaking in QED for any arbitrary weak interaction. The essence of this magnetic 
catalysis is that electrons are effectively 1+1 dimensional when their energy is much less 



than the Landau gap \j\eB\ what was pointed out recently in Refs.[^ ||. The lowest Landau 
level (LLL) plays here the role similar to that of the Fermi surface in the BCS theory of 
superconductivity, leading to dimensional reduction in dynamics of fermion pairing. 
The dynamical mass of fermions (energy gap in the fermion spectrum) is: 



rridyn — CVeB exp 



a 



where the constant C is of order one and a = e'^/An is the renormalized coupling constant. 

The effect of magnetic catalysis was studied also in Nambu-Jona-Lasino (NJL) models 
in 2+1 and 3+1 dimensions 0, it was extended to the case of external non-abelian 
chromomagnetic fields and finite temperatures H, ||], curved spacetime DTOf, as well as to the 



super symmetric NJL model [|TT|], confirming the universality of the phenomenon. 



We emphasize that we will work in the conventional, weak coupling, phase of QED. That 
is, the bare coupling a^^\ relating to the scale fi = A, where A is an ultraviolet cutoff, is 
assumed to be small, a^^^ -C 1. Then, because of infrared freedom in QED, interactions 
in the theory are weak at all scales and, as a result, the treatment of the nonperturbative 



dynamics is reliable. 

An important question of the chiral symmetry restoration in QED at finite temperature 



was addressed in the recent work of Lee, Leung and Ng who have obtained for the critical 



temperature ~ ■^y2Tr\eB\. However, their Tc can be considered only as a rough upper 
estimate. We will show that the correct estimate for the critical temperature is Tc ~ mdyn 
with TUdyn given by (|lD. 

The Lagrangian density of massless QED in a magnetic field is: 

1 Ft 



C = — F^T.^ + 
4 ^ 2 



^,i^rD,)^|J , (2) 



where the covariant derivative is 



D, = d,~ ie{Al^' + A,), = (o, |a:i, o) . (3) 

The Lagrangian density (|^) is chiral invariant ( we will not discuss the dynamics related 
to the anomalous current js^, in any case it is not manifested in quenched approximation 
dealt with in present consideration). When the chiral symmetry is broken there appears a 
gapless NG boson composed of fermion and antifermion. The dynamical mass for a fermion 
can be defined by considering the Bethe-Salpeter (BS) equation for NG boson [||, or the 
Schwinger-Dyson (SD) equation for the dynamical mass function [|13|, We consider here 
the homogeneous BS equation for NG bound state which takes the form: 

Xafi{x,y;P) = j d^Xid%id'^X2d%2GaaAx,Xi)Ka^f3^;a2(i2i^iyi,X2y2) 

■ Xa2i32ix2,y2;P)Gf3,p{yi,y), (4) 

where the BS wave function Xap{,x,y; P) = {0\TilJa{x)ipi3{y)\P) , and the fermion propagator 
Gai3{x,y) = {0\TilJa{x)ilj/3{y)\0) . Note that though the external field A^^^ (H) breaks the 
conventional translation invariance, the total momentum P is a good, conserved, quantum 



number for neutral channels |T^, in particular, for this NG boson. Since, as will be shown 



below, the NG boson is formed in the infrared region, where the QED coupling is weak, one 



can use the BS kernel in leading order in a : 

^«i/3i;a2/32(a;il/i,a;2?/2) = -47ria7^^„27^2/3i^^A*<^(?/2 - a;2)5(xi - X2)5(yi - ^2), (5) 
where the photon propagator 



(A is a gauge parameter). Then the BS equation takes the form: 



(6) 



S{x, xi)-f^x{xi,yi] P)YS{yi,y) 



l^f^uiyi -xi), (7) 



where, since we are working to the lowest order in a, the full fermion propagator G{x, y) is 
replaced by the propagator S* of a free fermion (with the mass m = rridyn) in a magnetic field 



S{x, y) = exp ( |(x - + y) ) S{x - y), 



where the Fourier transform of S is 
S{k) 



00 




J (is exp 


is ^ 








(fcV - + mdyn)il + 7'7'tan(eSs)) - k±-f±il + tan\eBs)) 



.2 .2 ,2 tan(egs) , ' 

,1.2 



(9) 



Here k± = (^1,^2)5 7± = (71572)- Using the variables, the center-of-mass coordinate, R 
{x + y)/2, and the relative coordinate, r = x — y, equation (^ can be rewritten as 

r — ri' 



Xa(iiR,r;P) = -Ana / (fRi(fri 
'r — Ti 



S[R-Ri + 



s 



R + Ri 



I)^,(-ri)exp -ie{r + nrAf{R-R^] 



(10) 



Here the function XapiR, P) is defined from the equation 

Xap{x,y;P) = {0\TMx)My)\P) = exp[zerM^^*(i?)]xa/3(i?, r; P). 



:ii) 



It is important that the effect of translation symmetry breaking by the external field is 
factorized in the phase factor in Eq.(|ll]) and equation ( p!OD admits a translation invariant 



solution, r;P) = exp{—iPR)x{r;P). Henceforth we will consider the case P = 0, 
corresponding to NG bound state. Then, transforming this equation into momentum space, 
we get: 



-Ana 



exp(— 



^ J a/3 

where p\\ = ijP.p^), p± = {p^,P^)- 



The crucial point for the further analysis will be the assumption that ^ y k-^l 

and that the region mostly responsible for generating the mass is the infrared region with 



k <^ y\eB\. The assumption allows us to replace the propagator S in Eq.(|T^ by that 
projected into LLL. In order to show this, we recall that the energy spectrum of fermions 
with m = rrifiyn in a magnetic field is 



E„{ps) = ±\/m%^ + 2\eB\n + pl 



n = 0,1,2,... 



(13) 



(the Landau levels). The propagator S{p) can be decomposed over the Landau level poles 
as follows I, |16| : 



S{p) = i exp 



pI - pI - mly^ - 2\eB\n 



(14) 



with 



DnieB,p) = {fi^'-p'r + mdyn 



1 - z7'7'sign(eP))L„ 2 



pi 



l + z7Vsign(eP))L„_i 2 



Pj 
\eB\ 



eB\) 

+ 4(pY+pV)^Li (2- 



Pj 
eB\ 



where L„(x) are the generalized Laguerre polynomials (L„ = L"^(x) = 0). Eq.(|l^) 
implies that at p^ rridyn ^ the LLL with n = dominates and we can take 



S{p) 



I exp( 



\eB\ p\-mly, 



1 — i7^7^sign(ei?)). 



(15) 



where p\\ = p^'-f^ — p^'y^ and p'j = {p^Y ~ {P^Yi Eq.(]T^ transforms into the following 
one (further for concretness we assume eB > Q): 



dyn 



■ p{k\l,k^)'^^^Y{l-^^V)V,.{kll-Pll,k^-A^), (16) 
where 

p{p\\,P±) = {P\\ - mdyn)x{p){P\\ - mdyn), (17) 

and / = |ei?|^^/^ is the magnetic length. Eq.(p!6|) implies that pip\\,p±) = exp(—Pp\)Lp{p^^), 
where f{p\\) satisfies the equation: 

.(.„) - ^/A,(i-7V)r|±^.(^,,)|±^ 

• ^^J.(^||-P||), (18) 



J d^k^expi ^)V^,{k\i -p\\,k^). (19) 



Thus the BS equation has been reduced to a two-dimensional integral equation. Of course, 
this fact reflects the two-dimensional character of the dynamics of electrons from LLL, that 
can be explicitly read from Eq . (p^) . 

We emphasize that the dimensional reduction in a magnetic field does not affect the 
dynamics of the center of mass of neutral bound states (in particular, this NG boson). 
Indeed, the reduction 3 + 1— 1 + 1 in the fermion propagator, in the infrared region, reflects 
the fact that the motion of charged particles is restricted in directions perpendicular to the 
magnetic field. Since there is no such a restriction for the motion of the center of mass 
of neutral particles, their propagator must have a (3 + 1) -dimensional form. This fact was 
explicitly shown for neutral bound states in NJL model in a magnetic field, in 1/Nc expansion 



0], and for neutral excitations in nonrelativistic systems [|T^]. This in particular implies 
that, notwithstanding the dimensional reduction in a magnetic field, the phenomenon of 
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spontaneous chiral symmetry breaking in QED does not contradict to the Mermin-Wagner- 



Coleman theorem [IS] forbidding the spontaneous breakdown of continuous symmetries at 
D = l + l. 

Henceforth we will use Euclidean space with = —ik^. In order to define the matrix 
structure of the wave function (p{p\\) of the NG boson, note that, in a magnetic field, there 
is the symmetry 5*0(2) x 5*0(2) x P, where the 5*0(2) x 5*0(2) is connected with rotations 
in X1-X2 and X3-X4 planes and V is the inversion transformation X3 — X3 (under which 
a fermion field transforms as ^ ^75 73 '?/')• This symmetry implies that the function '^{p\\) 
takes the form: 

V'(Pll) = 75(^ + ^7172^ + P\\C + 27i72P||-D), (20) 

where = ^373 +P474 (with 7^ anti -Hermitian in Euclidean space) and A, 5, C and D are 
functions of p|. 

In Feynman gauge 

dx exp(— /^a;/2) 



+ X 



and, substituting expansion ( pO]) into equation (0), we find that B = —A, C = D = 0, i.e., 
ip{p\\) = ^75(1 — "^7172), and the function A satisfies the equation 

a f (fkA{k) y dx exp{—xP /2) 

J k^ + m2 J {k-pf + x ^^^^ 
^ 

(henceforth we will omit the symbol || from p and k). Eq.(^T]) coincides with the equation for 
the dynamical fermion mass function obtained with the help of SD equation for a fermion 
propagator in [jl2|, |14|. As was shown in (see Appendix C), in the case of weak coupling 
a, the function A{p) remains almost constant in the range of momenta < p'^<l/f' and 
decays like outside that region. To get an estimate for m^y„ at a << 1, we set the 
external momentum to be zero and notice that the main contribution in the integral on the 
right hand side of Eq.(^) is formed in the infrared region with /c^<l//^. The latter validates 
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in its turn the substitution A{k) — > A{0) in the integrand of (|21|), and we finally come to the 
following gap equation: 

A O ^ / , / rk—^^ 22 

27r"' J k'' + m%y^ J + x 

i.e., 

a f ax exp[—ax) , _ _ ?«dj;ri' 





The main contribution in (|23| ) comes from the region x<l/a, thus at a << 1 we get 

therefrom the expression (|l]) for m^j^^ follows. The exponential factor in nidyn displays the 
nonperturbative nature of this result. It can be shown also that the expression (|I]) for the 
dynamical mass is gauge invariant 

More accurate analysis 0, which takes into account the momentum dependence of the 
mass function, leads to the result 



rridyn ^ CJ\eB\ exp 



IT IT 

"2 V 2^ 



(25) 



Notice that the ratio of the powers of this exponent and that in Eq.(|T|) is 7r/2-\/2 ~ 1.1, thus 
the approximation used above is rather reliable. 

To study chiral symmetry breaking in an external field at nonzero temperature we use 
the imaginary time formalism. Now the analogue of the equation (^) (with the replacement 
fT^dyn m?{T) in the denominator) reads 

A( \--T \^ 1 dkA{uJn,k) J dxexp{-l'^x/2) 

A[ur,,,p) ^1 l^j^ u;l + k^ + m^T) / (^„ - u^^f + {k - pf + x' ^''^ 

where ujn = TiT{2n + 1) are Matsubara frequencies. 

If we now take n' = 0, p = in the left hand side of Eq.(p6D and put A{ujn, k) ^ A{ujq, 0) = 

const in the integrand, we come to the equation 

1 - 7 7 dxexp{-l^x/2) 

~ 2^ /, ,2 I 7.2 I ^2frr\ / /, , , , \2 I u2 I ^- y^^') 



n n^oo ujl + k"^ + rn?{T) J {ujn - uJqY + k"^ + x 



After evaluating the sum in (pTj), we are left with the equation 

dkdx exp[-fx/2] J {irTf + x- m'^{T) 



a 

71 







[(7rT)2 + x- m2(T)]2 + {2nTy{k^ + m^{T)) 



2T / ' yFT^ 

The equation for the critical temperature is obtained from 



2T 

putting m{Tc) = 0: 



(2^ 



rffcda;e-2^("^^^)' [ 1/4 + x 



1 = - i 

' 



tanh (vrA;) H — ^ coth (tta/Zc^ + 



where we also switched to dimensionless variables x (27rTc)^x and A; 2TxTck. 



(29) 



By assuming smallness of the critical temperature in comparison with the scale put by 
the magnetic field, T^^ -C 1, we see that the double logarithmic in field contribution in 
Eq.(p9D comes from the region < x<l/2(7rTc£)^, 1/Tr<k < 00. Simple estimate gives: 

l/2(7rT,^)2 00 



a 

1 ~ - 

TT 



dx 



dk 







I/tt 



1/4 + x 1/4 -X 

[1/4 + x]2 + F [ k ^ yFT^_ 



l/2(7rT,^)2 



TT 



dx 



1 



2(1/4 + x; 



■log(l + (l/4 + x)27r2) + 



1/4 -X 



(1/4 + X + |l/4 - x|) Jl/7r2 + (1/4 + x)2 

• log , 

(1/4 + x)y^l/7r2 + x + |l/4 - xl/vr 

Thus, for the critical temperature, we obtain the estimate: 



a 
An 



(l/4 + x)|l/4-x| 
1 



log^ 



2(7rTe^ 



(30) 



eB I exp 



mdyn{T = 0), 



(31) 



where rridyn is given by (P. The relationship Tc ~ rridyn between the critical temperature and 
the zero temperature fermion mass was obtained also in NJL model in (2+1)- and (3+1)- 
dimensions [§, ^. The constant C, in the relation = Crridyn, is of order one and can be 
calculated numerically. We note the photon thermal mass, which is of the order of ^/aT 
p!9| , cannot change our result for the critical temperature. As is easy to check, the only 
effect of taking it into account will be the shift in x for a constant of the order of a in the 
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integrand of (|30|). However, such a shift is absolutely irrelevant for our estimate (|3TD. On 
the other hand, there are relevant one-loop contributions in the photon propagator due to a 
magnetic field. Taking them into account (that corresponds to the so-called improved ladder 
approximation) , we get an expression for rrifiyn of the form (|l|) with the replacement a —>■ a/2 

i- 

In conclusion, let us discuss possible applications of this effect. One potential application 
is the interpretation of the results of the GSI heavy-ion scattering experiments in which 
narrow peaks are seen in the energy spectra of emitted e~^e~ pairs. One proposed explanation 



1 21 1 is that a very strong electromagnetic field, created by the heavy ions, induces a phase 
transition in QED to a phase with spontaneous chiral symmetry breaking. The observed 
peaks are due to the decay of positronium-like states in this phase. The catalysis of chiral 
symmetry breaking by a magnetic field in QED, discussed in this report, can serve as a 
toy example of such a phenomenon. In order to get a more realistic model, it would be 
interesting to extend this analysis to non-constant background fields 

An interesting application of the magnetic catalysis in astrophysics, such as the cooling 
process of neutron stars due to enhanced Primakoff process in high magnetic field, was 
mentioned recently in p2| . 

Yet another application of the effect is connected with the role of chromomagnetic back- 
grounds imitating gluon condensate in the QCD vacuum H, 0, . 
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